Multi-interaction mean-field renormalization group 



scC. N. Likos^^^ and A. MaritarP'^ 

Dipartimento di Fisica Teorica, Universita di Trieste, 
Strada Costiera 11, 1-34014 Grignano (TS), Italy 

Istituto Nazionale di Fisica della Materia (INFM), 
Scuola Internazionale Superiore di Studi Avanzati 
and INFN, Sezione di Trieste, via Beirut 2-4, 1-34013 Trieste, Italy 



ABSTRACT: We present an extension of the previously proposed mean-field renor- 
malization method to model Hamiltonians which are characterized by more than just one 
type of interaction. The method rests on scaling assumptions about the magnetization of 
different sublattices of the given lattice and it generates as many flow equations as coupling 
constants without arbitrary truncations on the renormalized Hamiltonian. We obtain good 
results for the test case of Ising systems with an additional second-neighbor coupling in 
two and three dimensions. An application of the method is also done to a morphological 
model of interacting surfaces introduced recenlty by Likos, Mecke and Wagner [J. Chem. 
Phys. 102, 9350 (1995)]. 

PACS Nos.: 64.60.Ak, 64.60.Fr, 05.70.Jk 



Submitted to: Physical Review E 



I. INTRODUCTION 



One of the simplest techniques for deriving renormahzation group recursions is the so- 
called mean-field renormahzation group (MFRG) of Indekeu et al} The MFRG has proven 
to be an efficient method for the calculation of critical couplings and exponents of a wide 
variety of lattice models. The range of applications of the MFRG extends beyond the study 
of bulk critical phenomena in classical systems, to include quantum models and surface 
criticalities; moroever, extensions of the MFRG ideas to dynamical critical phenomena 
have also been presented. ^'"^ 

The most appealing features of the MFRG arc its simplicity and its direct connection 
with the classical mean-field ideas; whereas mean-field theory is a simple approach to the 
study of the phase behavior of various models in statistical mechanics, its most serious 
drawback is that it predicts the wrong (classical) critical exponents. The MFRG offers a 
way to remedy this deffciency while maintaining some basic notions from the mean-ffeld 
approximation (e.g. the effective field.) 

We shall now summarize the main ideas of the original MFRG. To illustrate the 
method, let us consider the nearest-neighbor Ising model which is described by the Hamil- 
tonian: 



where the second sum is carried over nearest-neighbor sites on a given lattice of dimension 
d. Consider now two clusters of interacting spins, containing N and N' sites, with N' < N. 
Suppose that the surrounding spins of these clusters is fixed at the values b (|6| < 1) for 
the large cluster and b' {\b'\ < 1) for the small one and evaluate the magnetizations per 
site for both clusters, m and m' respectively. Clearly, m and m' depend on the coupling 
constant J, the magnetic field h and the boundary values b and b'. Setting m( J, h,b) — b or 
m'(J, h, b') = b' we would obtain two different types of mean-field equations, which could 




(1.1) 



<ij> 
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then be solved self-consistently and the criticahty would be identified with the 'bifurcation 
point' of any of the two equations for = (each of the two equations yields a difi"erent 
approximate value Jc for the critical coupling, of course.) However, this approach leads to 
classical critical behavior. Instead, in the MFRG one uses the two cluster magnetizations 
to define a mapping (J, /i, h) — > ( J', /i', h') by requiring 



together with 



m'(j', h\ h') = s'^-y^miJ, h, b), 



(1.2) 



(1.3) 



to hold to leading orders in h and h. In Eqs. (1.2) and (1.3) above, the rescaling factor s 
is defined as 

s = {N/N')^'^. (1.4) 

Setting h = h' = 0, expanding both sides of (1.2) to linear order in b and b' and using (1.3) 
we obtain the RG flows for the coupling constant in the form 

dm'{J',h',b') dm{J,h,b) 



h=0,b=0 



(1.5) 



h'=o,b'=0 db 

Once the fixed points J* of the iteration have been found from (1.5) the thermal and 
magnetic exponents yt and yj^ are determined through the relations 

InAf 



yt = 



Ins ' 



where 



and 



'd'^m{J, h, b)\ fd^m'iJ', h', b')\-i 



dbdJ 

dm'{J',h',b') 



Ko^m [J , h ,b )\ 
WdJ' / 



J=J'=J^ ,h=h'=0,b=b'=0 



dh' 



_„d-2yh 9m{J,h,b) 
~' dh 



J=J^,h=0,b=0 



(1.6) 



(1.7) 



(1.8) 



J'=J^,h'=0,b'=0 

As a rule, the MFRG yields better estimates for the critical coupling Jc = J* than it 
does for the critical exponents. However, the latter can also be significantly improved in 
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a number of ways. Slotte^ has shown that a redefinition of the rescahng factor s leads to 
better values for the thermal exponent than the 'naive' definition of s, cq. (1.4) above. On 
the other hand, Indekeu et al.^ have succeeded in improving the MFRG critical exponents 
in a more systematic way, by considering three clusters and introducing an additional 
surface critical exponent y/,^; this approach has led to a unified approach to bulk, surface 
and corner critical behavior (see also Ref. 2.) 

The remaining of the paper is organized as follows: in Section II we motivate and 
present a new, extended version of the MFRG. In Section III we test the new approach by 
applying it to the square, simple cubic and body-centered cubic Ising models with crossing 
bonds, obtaining very good results for the fixed points and the overall phase behavior of 
these models. In Section IV we then apply this RG scheme to a morphological Hamiltonian 
of interacting surfaces.^ In Section V we summarize and conclude. Since the derivation of 
the RG fiows for the morphological model is a bit lengthy, we outline the main steps in 
the Appendix. 



II. EXTENDED MFRG 



Let us consider, to begin with, what happens if we want to renormalize, by means of the 
simple MFRG-scheme, a Hamiltonian in which the even interaction part {i.e. excluding the 
magnetic field and having a remainder which is invariant under — > — s^) contains more 
than just one type of interaction and thus more than one coupling constant. Denoting by 
J = (Ji, J2, . . . , Jp) the z^-dimensional coupling constant of the Hamiltonian, the MFRG 
which is based on the scaling of a single order parameter (the bulk magnetization) always 
leads to a single flow equation, namely one of the form (1.5) with J and J' replaced by 
their multidimensional counterparts J and J' respectively. Clearly, in order to specify the 
flows one needs as many equations as coupling constants in the Hamiltonian, so this single 
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equation is not sufficient. In the vast majority of MFRG applications to many-couplings 
Hamiltonians, this equation was used in the following sense: a 'fixed-point requirement' of 
the form J = J' was made, which led to an equation of the form /(J) = 0, identified as 
the equation which defines the critical surface of the model. Although such an approach 
has been widely used,^ and gives qualitatively correct results, it has the obvious drawback 
that it yields infinitely many 'fixed points' (each point on the critical surface is 'fixed') 
and thus detailed information on the critical exponents is lost in this scheme. Moreover, 
the approach only serves to define the critical surface alone and not to determine the 
way in which each of the individual couplings Jq,, a = 1,2, ... ,1/ flows in the presence 
of the others. Therefore, such an approach is useful only in the sense that it provides us 
with better estimates of the critical surface of a given model than the simple mean-field 
approximation, requiring roughly the same amount of computational effort; but it does 
not offer a means for the renormalization of the Hamiltonian. 

There have been a few attempts to go beyond the above limitations of the MFRG: de 
Oliveira and Sa Baretto^ studied the two-parameter Hamiltonian of the Ashkin- Teller (AT) 
model^ and identified two order parameters in the problem. Making scaling assumptions 
about the order parameters of small and large clusters, they were then able to obtain fiow 
diagrams of the usual type, finding isolated fixed points, critical exponents etc. However, 
the AT model is one in which two different types of spin variables live on the lattice sites, 
and this makes the introduction of two types of 'magnetization' natural in this model. An 
approach which has some similarities to ours was presented by Plascak^ for the case of the 
two-dimensional Ising model with crossing bonds. This approach is based on the selection 
of appropriate pairs of order parameters, depending on the region of the parameter space 
considered. However, this method leads to some ambiguities in the fiows along the axis of 
vanishing first-neighbor coupling.^ On the other hand, as we show below, our approach is 
free of such ambiguities. 
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We first present the basic conditions to be satisfied for the implementation of the 
new method, at present. The types of models on which our approach is applicable must 
satisfy the following requirements (the discussion of more general cases is postponed for 
the concluding Section): 

(i) the interaction part of the Hamiltonian must be even, and it must contain only two 
coupling constants, call them J and K\ 

(ii) the model must be defined on a bipartite lattice, and 

(iii) there must be a region of the parameter space for which there are only two kinds 
of ground states, a ferromagnetic (FM) and antiferromagnetic (AFM) one, separated by 
some borderline of stability. 

Once these prerequisites are satisfied, the approach proceeds as follows: consider, as 
in the original formulation, two clusters of N and N' sites with N' < N. (Hereafter, 
primed and unprimed quantities will refer to quantities pertaining to the small and large 
clusters respectively.) The clusters must be chosen in such a way that the two different 
sublattices of the given bipartite lattice are mutually equivalent in both clusters. Denoting 
the two sublattices by A and B and fixing the surrounding sublattice magnetizations to 
the values b'^ (bi) and 62 (^2) on the A and B sublattices respectively, we then obtain 
by the usual mean-field procedure expressions of the type m'j^i^^^{J' , K' ,h' ^h'-^^h'^ and 
mj^(^Q^{J, K, h,bi,b2) for the sublattice magnetizations of the small and big clusters. A 
mapping (J, K, h, bi, 62) — > ( J', K', h', b'^, 62) is now defined, in analogy with eqs. (1.2) and 
(1.3), through the requirements: 



m'^if, K', h', b[, b'2) = s^-y^fUAiJ, K, h, 61, 62), 



(2.1) 



and 



i= 1,2. 



(2.2) 
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Equation (2.2) simply expresses the intuitively appealing requirement that the boundary 
magnetizations scale the same way as the bulk ones, as in the original MFRG.-*^ We adopt 
the usual definition s = {N/N')^/^ for the rescaling factor. Setting again h = h' = 0, 
expanding both sides of (2.1) to linear order in b'-^ ^ h\^2 ^■nd using (2.2) we arrive at 
the flow equations for the coupling constants 

am^( J', K\ h', h') _ dniAi J, K, h, b 



h'=o,h'=0 dbi 



and 

dm'^{J', K', h', W) _ druAiJ, K, h, b) 



db'2 



(2.3) 



(2.4) 



/i'=0,b'=0 dhi 

where we have used b' and b as a shorthand for (6'^^, and (6i, h^). Due to the equiva- 
lence of the two sublattices, the individual subcluster magnetizations obey the symmetry 
mA{J, K, h, bi, 62) = mB{J, K, h, 62, 61) and the same for the primed quantities, thus the 
flows obtained if we use the B-sublattice magnetizations to perform the mapping are iden- 
tical to (2.3) and (2.4) above. Hence, the method generates exactly two independent flow 
equations, which are necessary and sufficient for the renormalization of the two-parameter 
Hamiltionian. 

After the fixed points J* = (J*,i^*) of the fiows have been found, the nonmagnetic 
eigenvalues Ai^2 ^-^d the critical exponents yi^2 = lnAi^2/lns, as well as the associated 
eigenvectors, are evaluated by the usual procedure of linearizing around J*. The remaining 
magnetic exponent is calculated through the relation 



dm'^{3',h',h') _ d-2y,dmAi3JiM 

— s 



dh' 



(2.5) 



j'=j^,/i'=o,b'=o dh 
Once more, due to the symmetry between and mentioned above, it is irrelevant 
which of the two is used for the evaluation of For the same reason, we can also work 
with the total magnetizations m' = m'^ + m'^ and m = rriA -\-'niBi derive yj^ from the 
scaling of the susceptibilty x = dm/dh through 

X' = s'^-'^y^x. (2.6) 



which will be useful later. We also note that, although the feature of working in the 
neighborhood of b' = b = is appropriate to the study of second-order transitions, the 
flow equations presented above are capable of producing also low-temperature flxed points 
with the associated y/j = cZ magnetic exponent which is the signature of a first-order phase 
change. 

Having developed the theoretical framework for our approximate renormalization tech- 
nique, we proceed in the following two Sections with specific applications. 



III. ISING MODEL WITH CROSSING BONDS 

In this section we apply the extended MFRG (EMFRG) to the case of the Ising model 
with first- and second- neighbor interactions (Ising model with crossing bonds.) The Hamil- 
tonian reads as 

i <ij> 

where the second sum is carried over nearest sites and the third over next nearest sites. 
The model satisfies requirement (i) of Section II; in order to meet the requirement of a 
bipartite lattice, we must restrict our choices: in two dimensions, we will consider the 
square (sq) lattice, which can be split into two square sublattices. Moreover, the model 
on the square lattice has three types of ground states (for h = 0): a doubly degenerate 
ferromagnetic (FM) one in the region {J > 0;K > —J/2}, a doubly degenerate anti- 
ferromagnetic (AFM) one in the region {J<0;K> J/2} and a four-fold degenerate 
one with super-antiferromagnetic (SAF) order in the remaining region K < —\J\/2. The 
SAF-ground states consist of alternating rows or columns of up and down spins. In order 
to satisfy requirement (iii) of Section II, we shall only consider the flows in the subspace 
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K >0 in which we have only the FM and AFM ground states seperated by the borderhne 
of stabihty J = 0. 

In three dimensions, we will study the simple cubic (sc) lattice which separates into 
two face-centered cubic (fee) sublattices, and the body-centered cubic one (6cc) which 
separates into two simple cubic sublattices. Again, for reasons similar to those of the two- 
dimensional case, we will only examine the flows in the subspace K > 0; for both the sc- 
and 6cc-models, the ground states are again FM for J > and AFM for J < 0. We will 
present the calculation in some detail for the sg^model only, and just describe the results 
for the sc- and 6cc-cases, since the essential characteristics of the flows are quite similar 
for all three lattices. 

A. Square lattice. For the sg^model, we take as the small cluster the N' = 2 spins joined 
by an elementary bond, and as the big cluster the A'" = 4 spins around an elementary 
plaquette of the lattice (fig. 1). These are the two smallest clusters we are allowed to 
consider which satisfy the requirement of equivalence of the sublattice-clusters. We take the 
points A,C,... forming the A-sublattice and the points B, D, . . . forming the -B-sublattice. 
With the boundary magnetizations fixed to the values b'^ (bi) and b'2 (62) for the A- and 
5-sublattices, we obtain the effective Hamiltonians: 

H'{SA, sb) = J'[SASB + KsAb2 + ss^l)] + 4K'{sAbi + 55^2) + h'{sA + sb), (3.2) 
for the N' = 2-cluster and 



H{sA, SB, sc, sd) = J SASB + sbsq + scsd + sdsa + 2[(s^ + sc)^2 + {sB + 
K 



+ 



SASC + SBSD + 3[(sA + sc)bi + {sb + S£>)&2] + h{sA + sb + sc + sq) (3.3) 

for the N = 4-cluster. Using eqs. (3.2) and (3.3) above, we find the sublattice magnetiza- 
tions per site as 

, _ sinh[(3 J' + AK') (6'i + ^'2) + 2/i'] - fi-^-^' sinh[(3 J' - AK') {b'^ -b'^)] 
™^ ~ cosh[(3J' + AK'){b\ + b'^) + 2h'] + cosh[(3J' - AK'){b\ - b'^)] 
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and 



AJ+2K 



sinh 



(4 J + &K){hi + 62) + 4/i +2 sinh(4J62 + QKhi + 2h)- 



-AJ+2K 



sinh 



(4J-6K)(6l-62)]}x 



(4 J + m){hi + 62) + 4/i +2 cosh(4J62 + + 2/i) + 



(4J-6K)(6l-62)] +2e-2^} \ (3.5) 



2 cosh(4 J61 + 6K62 + 2/i) + cosh 

Combining (2.3) and (2.4) with (3.4) and (3.5) above we obtain the (J, K) ( J', K') flows 
in the form 

3 J' + AK' - e-2-^'(3J' - AK') AJ + 6K + 12Ke-4^-2K _ e-8J(4j _ g^) 



1 + e-^J' 



1 + Ae-'^J-^K + e-8J + 2e-4./-4if 



(3.6) 



and 



3J' + 4K' + e-^-^'{3J' - 4K') 4J + 6K + SJe'^-^-^^ + e'^-^ {4J - 6K) 



(3.7) 



1 + 1 + 4e-4-/-2^ + + 2e-4-/-4^ " 

In the K > 0-subspace, the flows (3.6) and (3.7) have the fixed points {J^,K^) presented 
below. 

(i) A stable low-temperature fixed point Li = (+00, +00) corresponding to the FM ground 
state. 

(ii) A stable low-temperature fixed point L2 — (—00, +00) corresponding to the AFM 
ground state. 

(iii) A mixed low-temperature fixed point L3 = (0, -|-oo) which attracts along the J = 
direction only and repels in all other directions. This fixed point corresponds to the four 
degenerate FM and AFM configurations that are the ground states of the model which, 
for J = 0, reduces to two decoupled seising models with nearest-neighbor coupling K. 

(iv) A critical point Ci = (0.1590,0.1499) which represents the ferromagnetic Ising criti- 
cality. 
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(v) A critical point C2 = (—0.1590,0.1499) which represents the antiferromagnetic Ising 
criticahty. 

(vi) A multicritical point C3 = (0, 0.3465) which represents the criticahty of the uncoupled 
sg^Ising models described above, and 

(vii) the stable, high-temperature fixed point P — (0, 0) representing the paramagnetic 
phase. 

The flows in the K >0 subspace are shown in fig. 2. We note that some of the fiows 
emerge from two unstable fixed points located at J = —2K, J — > +00 and J = 2K, J — > 
—00, located at the borderlines of stability of the ground states, which are not shown. As 
a first remark, we note that the flows are very similar to those obtained in the past by the 
use of other RG techniques, ^^'^^ displaying reflection symmetry about the J — axis and 
a cusp of the critical surface at the point C3. However, the RG flows are obtained here in 
a much simpler way which is also readily generalizable to three-dimensional models in a 
very straigthforward fashion. The flows are also similar to those obtained in earlier work 
by Plascak,^ but with the signiflcant difference that here we avoid the ambiguities in the 
flows along the J = axis found in that work. In particular, we are capable to obtain 
the multicritical flxed point C3 on the J — axis which was missing in that approach 
and on which, as we demonstrate below, certain important relations between the critical 
exponents are satisfied to a reasonable degree of accuracy. On the other hand, an obvious 
deficiency of our approach in its present formulation is that the theory is unable to locate 
the SAF-critical fixed point at J = 0, < since the way of separating the sublattices 
makes it inapplicable at present to the region in which the SAF-fixed point lies. 

For the calculation of the magnetic exponent ?//, , a special treatment is required for 
the AFM fixed points. As can be seen from Eq. (2.6), when h is relevant, i.e. y/j > 
(and > d/2 of course) we have x' < X- the other hand, when h is irrelevant and 
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thus y/j < we obtain the inequahty x' > %. However, the last inequahty itself only 
implies y/, < d/2, and not necessarily yj^ < 0. Indeed, the sign of y/j is very sensitive to 
the choice of the rescaling factor s. In order to avoid complications related to the choice of 
this factor, and keep the discussion as simple as possible at present, we will characterize 
the magnetic field as irrelevant whenever we obtain x' > X-i ^^'^ this is indeed what we get 
at the AFM-criticality fixed point C2. 

For the low-temperature fixed points Li and L3 we obtain the correct magnetic ex- 
ponent = d, a manifestation of a first-order phase transition. For the point L2 the 
magnetic field is irrelevant, whereas at the paramagnetic fixed point P we find correctly 
the result yj^ — d/2. The locations of the T > critical points and the corresponding 
critical exponents are summarized in Table I. 

The critical lines intersect the K = Q axis at the points Jc = ±0.336, which is the 
estimate from our theory for the FM and AFM critical couplings of the nearest-neighbor 
Ising model, to be compared with the exact result"'^'^ Jc = ±0.441. This result is of the 
same quantitative accuracy with that obtained from the original MFRG with a 2 ^ 1 
mapping, namely Jc = 0.346; the latter is also the value of = Kc at the fixed point C3 
which represents the FM criticality of the uncoupled sq^Ising models with nearest neighbor 
interactions only. The thermal exponent yi = 0.760 (at points Ci and C2) is in reasonable 
agreement with the exact result yt = 1, and the magnetic exponent y/j = 1.487 differs from 
the exact one y/^ = 15/8 by about 20%. 

At the multicritical fixed point C3 we have an exponent yi = 1.073 with the cor- 
responding eigenvector along the (1, 0)-direction and an exponent 1/2 — 0.600 with the 
corresponding eigenvector along the (0, l)-direction. This allows us to identify the latter 
with the thermal exponent yt of the uncoupled Ising lattices, and we call yi = ys to comply 
with the terminology introduced in Ref. 12. As a first remark, we note that yt differs from 
the thermal exponent yi for the points Ci and C2 but this is not surprising in view of the 
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fact that our method is a 4 — > 2 mapping in general, but for J = it reduces effectively to 
a 2 ^ 1 mapping; different degrees of approximation in the mappings yield different esti- 
mates for the critical exponents, as expected. On the other hand, the magnetic exponent 
Ufi = 1.415 at C3 is not too different form = 1.487 at Ci and C2. According to van 

Leeuwen^^ the exponents ys and yt must satisfy the relations: 

ys = 2yh-d (3.8) 

and 

Vs/yt = Tising- (3-9) 

Using the results in Table I, we find that the rhs of (3.8) is equal to 0.830, whereas the Ihs is 
1.073, which shows once more that the critical exponents are not evaluated very accurately 
in the MFRG, at least with small clusters and with the naive definition of s; however, the 
ratio ys/yt is equal to 1.789, and the exact result for the rhs of (3.9) is 7ising = 1.75, which 
shows that ratios of the critical exponents come out rather accurately in this approach, 
since they are independent of the precise definition of the rescaling factor. 

We now extend the approach to three dimensions, studying the flows for the simple 
cubic and body-centered cubic models, always in the subspace K >0. 

B. Simple cubic lattice. For the sc-lattice the two clusters are again the elementary 
bond {N' = 2) and the elementary plaquette {N = 4), once more the smallest possible 
clusters which satisfy the requirements of the approach. The low- and high temperature 
fixed points Li, L2, and P are identical to the sq-case, also with the correct magnetic 
exponents. Once more, we find the two critical points Ci and C2 and the multicritical 
point C3 whose coordinates and critical exponents are summarized in Table II. The overall 
flow pattern is identical to the sg^case. 

The critical lines now intersect the K = 0-axis at the points Jc — ±0.192 to be 
compared with the 'exact' result"'^'^ Jc = ±0.222. At the critical points Ci and C2 we 
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find the thermal exponent yi = 0.506 which differs quite a bit from the 'exact' result^* 
yt — 1.587; the magnetic exponent y/j = 1.776 is somewhat closer to the value-'-'^ y/^ = 2.485. 
The value of the point C3 gives an estimate of the critical coupling of the fee nearest- 
neighbor Ising model, since for J = the model decouples into the two independent 
fee sublattices of the sc-lattice. We obtain Kc = — 0.091 to be compared with the 
estimate-*^"^ Kc = 0.102. We note in passing that the MFRG inherits from the underlying 
mean-field ideas the charactersistic that the values of the critical couplings improve, for the 
same degree of approximation, as the coordination number of the lattice increases. The 
value of the thermal exponent, however, does not: at the point C3, yt = y2 = 0.344 which 
is worse than the value at Ci and C2; ideally, all these estimates should converge to the 
exact 3D-Ising values as one considers larger and larger clusters. The magnetic exponent at 
C3, yfi = 1.688, is again not too different from its value at Ci. Regarding the relation (3.8) 
we obtain 0.443 for the Ihs and 0.377 for the rhs, whereas the ratio ys/yt = 1-288 is not 
too far away from the 'exact' value 7ising = 1-238 (Ref. 15) required from eq. (3.9). Again, 
the ratio of the critical exponents is very accurate, although the individual exponents are 
not. 

C. Body-eentered eubie lattiee. Here, the small cluster is chosen to be the nearest- 
neighbor bond {N' = 2). For the large cluster, there are two AT = 4 choices that satisfy 
the requirements: the ASC-D-tetrahedron and the ADCE parallelogram shown in fig. 3. 
Both mappings yield identical low- and high-temperature fixed points, with the correct 
magnetic exponents. Moreover, they yield the same multicritical point C3 = (0, 0.203) but 
they differ slightly in the nonmagnetic critical exponents at C3 as well as at the locations 
and exponents of the critical points Ci and C2. The flow patterns are again very similar to 
those presented above for the other lattices. Although the critical fixed points are different, 
the eritieal lines from the two mappings are almost identical. This demonstrates that the 
method is robust (insensitive) to the choice of the large cluster. 
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The results are summarized in Tables Ill(a) and Ill(b). The value = 0.203 at the 
fixed pint C3 is an estimate for the critical coupling of the sc-Ising model with nearest 
neighbor coupling, again not too far from the exact value, 0.222, and the previous sc 
estimate, 0.192. The critical lines now intersect the K = 0-axis at the points Jc = ±0.140 
for the tetrahedron bond mapping and Jc = ±0.139 for the parallelogram — > bond 
mapping which are the estimates for the 6cc-Ising critical coupling, to be compared with 
the best estimate-*^"^ Jc = ±0.157. A comparison of the critical exponents with the exact 
ones shows that the tetrahedron — > bond mapping yields better results not only in terms 
of their numerical values, but also in the sense that it predicts almost identical values for 
the thermal and magnetic exponents for the points Ci and C3, as required by universality. 
The relations (3.8) and (3.9) are satisfied within an error of at most 10%. 

The above discussion demonstrates that the RG-method we propose here yields very 
satisfactory results for some standard 'test'-models, and thus it is a plausible technique for 
the renormalization of more complicated Hamiltonians. 



IV. MORPHOLOGICAL MODEL 



We consider in this Section a phenomenological morphological Hamiltonian on a bcc- 
lattice, introduced recently in order to model the phase behavior of microemulsions.^ Here, 
we outline the basic ideas in the derivation of the model, and refer the reader to Ref. 6 
for details. Let us consider a three-dimensional Bravais lattice with periodic boundary 
conditions, having N sites and volume V, whose Wigner-Seitz (WS) unit cells contain 
either bulk water or oil. In addition, the system contains amphiphilic molecules which 
are supposed to form an incompressible membrane in the interface between oil and water, 
defining in this way a Gibbs dividing surface between the two bulk phases. After choosing 
an orientation for this interface, a collection of water cells fixes uniquely an interfacial 
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pattern, and vice versa. We denote a cell to be 'occupied' if it contains water, and 'empty' 
if it contains oil. The morphological features of each pattern can be described by the 
Minkowski functionals; in three dimensions, these are the geometric invariants: covered 
volume, surface area, integral mean curvature of the interface, and the combinatorial Euler 
characteristic X of the pattern. We have X = (number of disconnected components)- 
(number of handles)+ (number of cavities) in three dimensions. In the present model the 
interfaces have no holes. The family of the Minkowski functionals is characterized by 
a theorem which asserts that any real-valued, additive, motion-invariant and continuous 
functional defined on the collection of the 2"^ configurations is a linear combination of 
the Minkowski functionals. ■'^^'■'^^ In order to deal with the statistical morphology of the 
interfacial membrane, we therefore take the Hamiltonian to be of the generic form: 

3 

n=J2haVa, (4.1) 
a=0 

where ha are energy parameters and Va are, within proportionality factors, the dimension- 
less Minkowski functionals as follows: denoting by V, A, A4 and X the covered volume, 
exposed area, integral mean curvature and Euler characteristic of the pattern formed by 
the full cells, we have Vq = V/{8V2l^), Vi = ^/[3(4^/2 + 2)^2], V2 = M/{6Trl) and Vs = X. 
The various factors arise from the definition of the Minkowski functionals via Steiner's for- 
mula, and their values for the hcc WS polyhedron (see the Appendix and Table I of Ref. 6.) 
The length scale I is the edge length of the WS-unit cell of the &cc-lattice, / = a\/2/4, where 
a is the hcc lattice constant (see fig. 3.) 

Using the property of additivity, one can derive concise expressions for the Minkowski 
functionals in terms of the occupation numbers tt^ = 0, 1 (itj = if the cell is occupied, 
= 1 if it is empty) and in Ref. 6 the Hamiltonian was written down explicitly in that 
representation. Here, it will be more useful to write the Hamiltonian in terms of 'Ising 
spin' variables, = ±1 with = 1 denoting a full site, and Sj = —1 an empty one {i.e. 
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Sj = 1 — 2ttj.) Setting s = (si, S2, ■ ■ ■ , S]\i), the final expression reads as 



'^(s) = ^aVa{s) = ho Y + J XI 
a=0 



hi 
4 



N- 



%/3 



4a/3 + 2 



E 



SlSj 



3(4a/3 + 2) 



E 



+ 



S^Sj 



+ 



h2 
4 



2^ ^iSjSfcS; . (4.2) 

The primed and double-primed sums are carried over nearest and next-nearest neighbor 
bonds, respectively. The triple-primed sum runs over isosceles triangles, two of the sides 
of the triangles being first neighbor bonds and the third being a second neighbor bond. 
Such are the ABC, BCD, ABD, and ACD triangles in fig. 3, for instance. Finally, the 
four-primed sum runs over tetrahedra whose faces are isosceles triangles as above, e.g. the 
ASC-D-tetrahedron of fig. 3. Under the interchange — > —s^, the integral mean curvature 
is odd, whereas the exposed area and Euler characteristic are even. 

We are going to deal exclusively with the case h2 — (no spontaneous internal curva- 
ture.) For ho = h2 = 0, the ground states (GS's) of the model are the following: a doubly 
degenerate ferromagnetic GS (oil or water, OW) in the region 



V3 



-hi;hs > -hi }; 



(4.3) 



2^3+ 1 

a doubly degenerate antiferromagnetic GS ('plumber's nightmare', PN) in the region 
and a four- fold degenerate 'droplet '-phase in the region 

The 'droplet' phase is realized when one of the two sc-sublattices has ferromagnetic and 
the other antiferromagnetic order. The subspace h^ > is covered completely by the OW 
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and PN-ground states, separated by the borderline of stability hi = {2^/3 + l)h^/\^; thus, 
according to the general requirements laid out in Section II, we are going to consider the 
flows in the subspace >0 only.^^ 

Dropping the uninteresting spin-independent constants from the Hamiltonian, and 
deflning: 

h=-^; J=^^ ; K = ^, (4.6) 

2 24(4^3 + 2) 16 ' ^ ' 

we arrive at the expression 

I II nil 

H = -pn = /i ^ + (6^3 J - 2K) ^ SiSj + 4J^ SiSj + -P^ ^ SjSjSfcS;. (4.7) 

i 

For the choice — the model reduces to a conventional bcc Ising model with an addi- 
tional second-neighbor coupling, namely 

H = h^Si + SiSj + aL^ SiSj, (4.8) 
i 

where L = P^/3hl/{16^/3 + 8) and a = 2/(3-\/3) is the ratio of second to first neighbor 
coupling. 

We are now interested in the renormalization of Hamiltonian (4.2). When Hq = h2 = 
a real-space RG will generate flows with at least three parameters corresponding to 
the three even interactions in (4.7), whereas the initial morphological Hamiltonian (4.1) 
contains only two, hi and (or J and K in (4.6), (4.7)) in the even subspace. However, 
after coarse-graining the effective Hamiltonian must still be additive, motion-invariant 
and continuous as the initial one. Thus, it must be again a linear combination of the 
Minkowski functionals, i.e. of the form (4.1). Our extension of the MFRG allows us to 
impose this requirement quite naturally, by assuming that the renormalized Hamiltonian 
for the AT'-cluster always has the form (4.7). 

Clearly, the model satisfles all prerequisites of Section II, and we can proceed with the 
EMFRG approach. We choose for the small cluster the ^S-bond {N' = 2) and for the 
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large one the ^BCD-tetrahedron {N = 4) shown in fig. 3. The points A,D,... form the 
^-sublattice and the points B,C,. . . form the 5-sublattice. Some details on the derivation 
of the flows are given in the Appendix. The flow equations read as 

(42^3 J' - lAK') sinh(6v^J' - 2K') + 24J'cosh(6y3J' - 2K') _ 

cosh(6\/3J' - 2K') ~ 

^ |g(24x/3+8) J-T/C [^72^3 + 40) J - 22K] + 2e-^ {AQJ - 2K)- 
g(-24x/3+8)J+9i^[^72V3 - 40) J - 26i^]| x 
|g(24V3+8)J-7i^ ^ ^ 2e"^-^+^ + e(-24V3+8)'^+9i^ 

and 

{A2^/^J' - UK') cosh(6v^/ - 2K') + 24/ sinh(6\/3 - 2K') _ 

cosh(6\/3J' - 2K') 

^ |e(24\/3+8)J-7Jr[(72V3 + 40)J - 22K] + 26-^(72^3 J - 24K)+ 
g(-24x/3+8)J+9i^[^72V3 - 40) J - 26X]| x 
|g(24V3+8)J-7K ^ ^ 2e~^-^+^ + e(-24V3+8)J+9K 

The above flows have in the subspace K >0 the fixed points {J^, K^) listed below. 

(i) A stable low-temperature fixed point Li — (-|-oo, -|-oo) representing the ferromagnetic 
(OW) ground state. 

(ii) A stable low-temperature fixed point L2 = (—00, -|-oo) representing the antiferromag- 
netic (PN) ground state. 

(iii) A mixed low-temperature fixed point L3 = (J*, 3-\/3</*) with J* -\-oo. This point 
attracts in the direction K = 3\/3J (i.e. /i3 = \/3/ii/(2\/3 -|- 1)) and repels in the other 
directions. It corresponds to the four-fold degenerate ground states of the model when 
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the couplings hi and have the ratio given above, and in which case the first-neighbor 
couphng vanishes; then, the OW and PN phases are all degenerate at T = 0. 

(iv) A critical point Ci = (1.656 x 10~^, 5.306 x 10~^) which represents the Ising ferro- 
magnetic criticality. 

(v) A critical point C2 = (1.407 x 10~^, 9.858 x 10~^) which represents the Ising antifer- 
romagnetic criticality. 

(vi) A multicritical point C3 = (2.785 x 10^^, 1.447 x 10"^) and 

(vii) the high-temperature, stable fixed point P = (0, 0) representing the paramagnet. 

The flows are shown in fig. 4(a); due to the particular representation chosen to show 
the flows, the fixed points Li and L3 coincide in this figure. We obtain the correct magnetic 
exponent y/j = (i at Li and L3 and = d/2 at P. The subspace — \/3/ii/(2\/3 -|- 1) 
is an eigenspace of the flows, with the fixed point C3 lying in this subspace. Flows that 
start in this space remain in it, running towards the point L3 if they start above C3 or 
towards the point P if they start below C3. Thus, in the special case of a model with 
vanishing first-neigbor coupling, the flows maintain that property. The parameter space is 
thus separated into three basins of attraction (see fig. 4(a)). Flows starting in the region 
enclosed in the 'triangle' formed by the two critical lines and the axis = run towards 
the P-point. Flows starting outside the triangle and on the right of the line P — C3 — L3 
run to the OW fixed point Li, always remaining in the region < \/3hi/{2\/3 + 1); and 
flows that start outside the triangle but on the left of the line P — — run initially 
close to L3, but eventually they turn around to end up in the PN fixed point L2, staying 
always in the subspace > ^/3hi/{2\/3 + 1). 

The locations of the T > 0-fixed points and the associated critical exponents are 
summarized in Table IV. We note that the method predicts identical values for the thermal 
exponent yi for both ferromagnetic and antiferromagnetic Ising critical points, as should. 
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The value yi = 0.691 is comparable to the value yi = 0.653 obtained for the bcc Ising model 
with crossing bonds using the same kind of mapping (tetrahedron — > bond) discussed in 
Section III (see Table 111(a).) As further evidence for the consistency of the method, we 
note that the critical line of the Ci fixed point intersects the = 0-axis at the point 
(3h1 = 2.019 (see fig. 4(b)). Accordingly, the critical coupling of the Hamiltonian (4.8) (to 
which the model reduces in the case = 0) is predicted to be Lq = 0.098. On the other 
hand, the line corresponding to a = 2/(3-\/3) intersects the critical line of the ferromagnetic 
fixed point Ci of the bcc Hamiltonian studied in Section III at a point whose abscissa is 
equal to 0.105. The two estimates differ by less than 7%. 

In order to obtain more detailed information about the phase behavior of the model, we 
have to invoke the results from other techniques as well, for example from the simple mean- 
field approximation. The reason is that the proposed RG scheme has certain limitations 
due to the feature of always expanding the cluster magnetization around vanishingly small 
expectation values (s^) of the surrounding spins. Consequently, not all of the points lying 
on the 'critical lines' correspond to true criticalities. Let us consider, for instance, the 
mean-field expression for the free energy per site which reads as^ 



where m = (s^). Expanding the free energy about m = up to 0(m ) and dropping the 
uninteresting constants we obtain 



The requirement of vanishing coefficient of the quadratic term identifies the line of diverging 
susceptibility, x""*^ = 0, and reads as 



Pf{m) = ^(1 - m2) + ^(-1 + 4m2 - 3m'^) + 
4 8 




(4.11) 




(4.12) 



(3hi-2 
2 



(4.13) 
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However, not all the points defined by (4.13) correspond to true criticalities which happen 
only when the coefficient of the m'^-term is positive. This requirement yields a tricritical 
point {(3h^{ , (3h^{) = (22/9,2/9). In fig. 4(b) we plot the line (4.13) along with the fiows 
and indicate the tricritical point by a dot. We note that below the dot the attractive line of 
the fixed point Ci almost coincides with (4.13). The true tricrical point of the model must 
occur at a temperature lower than the mean-field prediction, of course, but taking for now 
the estimate for tricriticality from the mean-field approximation for granted, we can give 
to the attractive line of the point Ci the following interpretation: the segment of the line 
below the dot corresponds to true criticality. However, the part of the solid line above the 
dot does not represent a line of critical points, but rather a line of diverging susceptibility 
at m = which, however, does not correspond to true criticalities because these are 
preempted the phase coexistence between oil- rich and water-rich phases. -"^^ The inability 
of the present RG scheme to distinguish between a critical line and a line of preempted 
criticalities lies in the expansions of the magnetizations about vansihingly small values of 
the boundary spins, i.e. in the assumption that the system is translationally invariant 
and has a vanishingly small bulk magnetization. Similar considerations also apply to the 
C2-'critcal line'. 

Another feature of the model found in the mean-field approximation and confirmed by 
extensive Monte-Carlo simulations^ is the possibilty of three phase coexistence between oil- 
rich, water-rich and a middle disordered phase, characterized as a microemulsion in view 
of its morphological (negative Euler characteristic) and structural (a peak of the struc- 
ture factor at nonzero wavevector) properties. The three phase coexistence is caused by 
the competition between a positive surface tension (3hi ~ 1 of the incompressible surfac- 
tant film, which tends to minimize the exposed area, and a positive topological potential 
(5h^ ~ 1 which encourages structures with a large surface area in order to accomodate 
many handles, and thus it assists the entropic tendency to disperse the amphiphiles. The 
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microemulsion has a finite correlation length ^ of the order of a few lattice constants, i.e. it 
displays short-range order (unlike the completely random mixture or 'paramagnet' in the 
magnetic language.) Clearly, since ^ is neither vanishing nor diverging, the microemulsion 
is not represented by a fixed point in the flow diagram. We must, therefore, resort once 
more to external information in order to identify the region in the flow diagram which 
represents stable microemulsion phases. According to earlier work,^ the middle phase 
occurs in a range of temperatures O.SO^kBT/h^ < 9/2, the upper limit being the tri- 
critical temperature, above which three-phase coexistence ceases to exist. On the other 
hand, below the lower limit the middle phase again disappears because it is replaced by a 
slightly disordered 'plumber's nightmare' phase. We take the combinations of the energy 
coefficients that give a three-phase coexistence as 'pointers' that indicate the region on 
the flow diagram where the microemulsion is stable. In flg. 4(b) we mark a few of those 
points, which are once more obtained form the mean-fleld approximation; the latter has 
been found to be relatively accurate in its predictions for the triple points, when compared 
to the simulation results.^ Referring to this flgure, we can now assert that the neighbor- 
hood of the filled triangles below the PC3 line, to the left of the solid line and above the 
tricritical point roughly defines the domain of thermodynamic stability of the microemul- 
sion. Clearly, upon coarse-graining all points in the microemulsion regime flow towards 
the paramagnet. In order to delimit the region of stability of the microemulsion more 
accurately one could, for example, study the small wavenumber behavior of the Fourier 
transform of the correlation function within the RG scheme. However, such a calculation 
lies beyond the scope of this work. 



V. SUMMARY AND CONCLUSIONS 



We have presented a generalization of the mean-fleld renormalization group method 
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which provides a way of renormahzing Hamiltonians with two couphng constants in the 
even interaction part. The method builds on ideas which are similar to those of the original 
MFRG, but it goes beyond the limitation of a single flow equation by employing scaling 
assumptions about suitably chosen sublattice magnetizations. We applied the method to 
several model Hamiltonians, obtaining satisfactory results for the flow diagrams. An ob- 
vious challenge for the future is the development of these ideas even further, so that we 
will be able to deal with Hamiltonians involving more than two parameters. A possible 
way to achieve this goal will be through a separation of the lattice into more than just two 
suitably chosen sublattices; however, since the EMFRG always yields as many flow equa- 
tions as sublattices, it may be necessary in some cases to augment the original Hamiltonian 
by a suitably chosen number of interactions and coupling constants until the number of 
couplings matches the number of sublattices. Then, the flows would be obtained in this 
enlarged parameter space, and one could look at the flows in the original, restricted domain 
by taking the appropriate 'cuts' in Hamiltonian space. 

The morphological model of Section IV includes all additive geometrical invariants 
whose thermal averages are extensive. The manifest additivity of the Hamiltonian is a 
sufficient, but not necessary condition for the thermodynamic requirement of extensitivity 
of the internal energy. Therefore, one should not exclude a priori terms in the Hamiltonian 
which are not additive; in particular, a non-additive 'curvature-square' term, which is 
employed in most current models of microemulsions, is missing from the Hamiltonian. In 
the original paper, ^ it was argued that the model deals with length scales exceeding the 
persistence length^^ where the above contribution (also called 'bending energy') can 
be omitted because the scale-dependent bending rigidity has been renormalized away. A 
non-perturbative renormalization of a Hamiltonian including the surface area and bending 
energy terms only was presented recently our model is complementary to that of Ref. 21 
in that it includes the Euler characteristic term, but not the bending rigidity. It would be 
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desirable, therefore, to start with a Hamiltonian that includes the Minkowski functionals 
and the bending energy and proceed with its renormalization, in order to see the crossover 
from the rigidity-dominated regime (for lengths scales below ^k) to the regime above 
where the thermal fluctuations dominate over the rigidity and the membrane is crumpled. 
We plan to return to this problem in the future. 
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Mobility Programme of the Commission of the European Communities, Contract No ER- 
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Here we outline the steps for the derivation of the flow equations for the morphological 
Hamiltonian of Section IV, eq. (4.7). Consider first the small cluster AB, N' = 2. Fixing 
all the surrounding ^(S)-sublattice spins to the value &^(&2)) using (4.7) and figure 3, 
we arrive at the effective Hamiltonian of the N' = 2 cluster of the form: 



ACKNOWLEDGMENTS 



APPENDIX 



H\sa,sb) = Ii'S'q + (6^3/ - 2K')S[ + 4/^ + K' S'^ 



(A.l) 



where 



S'q = + SB, 



(A.2) 
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and 



S[ = Sy^SB + 7(5^62 + ssb'i), 
S'2 = Q{sAh'i + SBh'2) 



S4 = 6b'ib2[sASB + 3(sA^2 + ^B^'i)]- 



(A.3) 
(A.4) 

(A.5) 



Simliarly, the effective Hamiitonian for the ABCD-clustei {N = 4) has the form 



His A, SB, sc,sd) = hSo + (6^3 J - 2K)Si + AJS2 + KS^^, 



where 



S{) = SA + SB + sc-\- sjj, 



(A.6) 



(A.7) 



^1 = SASB + sbSd + sdSc + ^csa + 6[(sa + S£,)62 + (s^ + sc)hi] , (A.8) 



and 



'S'2 = SASD + SBSC + 5[(s^ + SD)hi + {sb + 5(7)62] 

S'4 = s^SBSc'S£) + SBSc{sA + + sasd{sb + sc)^2+ 

3(s^SB + SBSD + SD^C + sCSa)^1^2 + SASD^i + 55^(761+ 
136i62[(s^ + SD)h2 + (sB + sc)hi]. 



(A.9) 



(A.IO) 



Using the above expressions, and ignoring terms of order (6')^ and iP' which do not con- 
tribute anything to the flow equations, we find the sublattice magnetizations per site as 



tu'a = |e6v^-^'-2^' sinh [[(42^3 + 24)/ - IAK'] {b[ + h'^) 



+ 



^^•^'+2^'sinh[[(42v^ - 24)/ - lAK']{h'2 - b[] 
^6V^J'-2K' ^^gj^ r[(42v^ + 24)/ - UK'] {b[ + b'2) 

-6V3J'+2K' sinh[[(42\/3 - 24)/ - UK']{b'2 - 61)] | 
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X 

+ 
-1 



(A.ll) 



for the small cluster, and 



mA = |e(24V3+8)^-7i^ ginh [[(72v^ + 40) J - 22K] {hi + 62)] 

2e~^sinli[(72v^J - 2AK)h2 + (40 J - 2K)6i] + 
g(-24V3+8) J+9i^ sinh [[(72^3 - 40) J - 26K] (62 - 61)] 
,(24V3+8) J-7i^ ^^gj^ J[(72^ + 40) J - 22K] {hi + 62) 
2e~^ cosh [(72V^J - 24iC)62 + (40 J - 2K)h^ + 
2e-^ cosh [(72\/3 J - 24K)6i + (40 J - 2i^)62 + 

2e-^J+K ^ e(-24V3+8)J+9i^ ^^31^ J[(72V3 - 40) J - 26K](62 - ^l)] } 



+ 



X 



+ 



-1 



(A.12) 



for the large cluster (where we have set h' = h = 0.) Differentiation of eqs. (A. 11) and 
(A.12) with respect to 6^ h'2 and &i, 62 leads to the flow equations (4.9) and (4.10) of 
the main text. The zero-field susceptibilities needed for the calculation of the magnetic 
exponent read as 



,6\/3J*-2Js:* 



X* — 



and 



X* 



cosh(6V3J* - 2K^) 



j(24V3+8) J.-7if* _^ 4e-^* + 26-8-^*+^* + ^{-2'iV2,+S)J^+QK^ ' 



(A.13) 



(A.14) 
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TABLE I. The T > 0-fixed points of the RG flows for the sg^Ising model with crossing 
bonds presented in the text [eqs. (3.6) and (3.7)] along with the corresponding critical 
exponents. 



{J*,K^) yi y2 2y,^ - d 

Ci (0.1590,0.1499) 0.760 -1.023 0.974 

C2 (-0.1590,0.1499) 0.760 -1.023 -0.037 

C3 (0,0.3465) 1.073 0.600 0.830 

P (0,0) -1.170 -0.830 
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TABLE II. The T > 0-fixed points of the RG flows for the sc-Ising model with crossing 
bonds and the corresponding exponents. 



Ci (0.0629,0.0580) 0.506 -0.424 0.552 

C2 (-0.0629,0.0580) 0.506 -0.424 -0.056 

C3 (0,0.0912) 0.443 0.344 0.376 

P (0,0) -0.966 -0.377 
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TABLE Ill(a). The T > 0-fixed points of the RG flows for the 6cc-Ising model with 
crossing bonds and the corresponding exponents obtained from the tetrahedron bond 
mapping. 







yi 


y2 






(0.0781,0.0791) 


0.653 


-0.848 


0.724 




(-0.0781,0.0791) 


0.653 


-0.848 


-0.061 


Cs 


(0,0.2027) 


0.901 


0.650 


0.788 


P 


(0,0) 


-0.667 


-0.789 
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TABLE Ill(b). Same as in Table Ill(a), but using the parallelogram — > bond mapping. 







yi 


y2 


2y/i - d 




(0.0994, 0.0494) 


0.426 


-0.800 


0.466 




(-0.0994, 0.0494) 


0.426 


-0.800 


-0.026 




(0,0.2027) 


0.846 


0.650 


0.788 


p 


(0,0) 


-0.320 


-0.789 
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TABLE IV. The T > 0-fixed points of the RG flows for the morphological Hamiltonian 
and the corresponding exponents. 







) 


yi 


y2 


- d 


Ci 


(1.656 X 10" 


"2,5.306 X 10-2) 


0.691 


-0.740 


0.646 




(1.407 X 10" 


-2,9.858 X 10-2) 


0.691 


-0.610 


-0.029 


Cs 


(2.785 X 10" 


-2, 1.447 X 10-1) 


0.728 


0.564 


0.516 


P 


(0,0) 




-0.667 


-0.789 
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FIGURE CAPTIONS 



FIG. 1. The N' = 2 and N = 4 clusters used for the remormahzation of the seising 
Hamiltonian with crossing bonds. The fluctuating clusters are shown by solid lines. The 
filled dots form the A-sublattice and the open ones the S-sublattice. 

FIG. 2. RG flows of the sg^Ising Hamiltonian with crossing bonds. The critical line is 
denoted by solid lines. 

FIG. 3. A local arrangement of the fecc-lattice and the bcc Wigner-Seitz unit cell. 

FIG. 4. RG flows of the bcc morphological Hamiltonian; (a) throughout the subspace 
> and (b) in more detail, in the neighborhood of the critical points. The lines along 
which the critical points Ci and C2 attract are denoted solid. The dashed line is the 
locus of vanishing inverse susceptibility from the simple mean-field approximation. The 
mean-field tricritical point is indicated by the dot. The filled triangles denote values of the 
energy parameters for which the mean-field approximation yields three-phase coexistence 
between oil-rich, water-rich and a middle, random, bicontinuous microemulsion. 
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